Laser Ion Acceleration Scaling Laws Seen in Muti-Parametric PIC Simulations 
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The ion acceleration driven by a laser pulse at intensity I — 10^" — 10^^ W/cm^ x {^ra/Xf' from 
a double layer target is investigated with multi-parametric Particle-in-Cell (PIC) simulations. For 
targets with a wide range of thickness I and density rie, at given intensity the highest ion energy 
gain occurs at certain critical depth of the target a = n^l, which is proportional to the square root 
of intensity. In the case of thin targets and optimal laser pulse duration, the ion maximum energy 
scales as the square root of the laser pulse power. When the radiation pressure of the laser field 
becomes dominant, the ion maximum energy becomes proportional to the laser pulse energy. 

PACS numbers: 52.38.Kd, 52.65.Rr, 29.25.-t, 41.75.Lx 



(N 



I 



CZ3 

• I— I . 
Oh 



> 

OS 
00 



in 
o 

o 



X 



The laser-driven ion acceleration can be of benefit 
for many applications, e. g., hadron therapy in oncol- 
ogy fast ignition of thermonuclear fusion by protons 
[2'| , production of sources for positron emitting tomogra- 
phy (PET) 0, conversion of radioactive waste 0|, etc. 
The radiation pressure dominant (RPD) regime of the 
laser ion acceleration can be a basis for a laser-driven 
heavy ion collider and proton dump facility for neutrino 
oscillation studies 6]. All these applications can be- 
come possible owing to invention of high-intense lasers 
0, that now are capable to produce pulses with inten- 
sity 10^" - W^YJ/cm^. 

Irradiation of solid targets by lasers with power from 
terawatt (TW) to petawatt (PW) and intensity up to 
10^1 W/cm^ resulted in generation of proton beams with 
maximum energy £max from 1 MeV up to 55 MeV 
Typical proton beam has picosecond duration and is 
emitted in quasilaminar fashion |^; it can have very low 
transverse emittance |lOj . Experiments indicate that the 
proton energy increases with the laser intensity and it 
depends also on the target (foil) thickness and the target 
composition. It turns out that for given laser intensity 
there is an optimal target thickness, at which the proton 
energy is maximum [llj. Thus the important question, 
what is the ion maximum achievable energy at given in- 
tensity and how it scales with increasing intensity? 

In this Letter we investigate the scaling laws of the ion 
acceleration driven by a laser pulse at intensity / — 10^" — 
10^2 (W/cm^)(/i,m/A)^ from a double layer target with 
multi-parametric Particle-in-Cell (PIC) simulations. The 
double layer target was suggested in Refs. [ij for a laser- 
driven high-quality ion beam generation that can be used 
in hadron therapy; this scheme was also demonstrated 
with three-dimensional PIC simulations in Refs. fl^ . 

In publications on the multi-MeV ion acceleration 
driven by an intense laser from foil targets, several 
regimes of acceleration can be found. For example, ions 
can be accelerated by varying electrostatic potential of 
the electron cloud expanding into vacuum, as shown in 
Refs. 01 with different models. In the case of thin foils, 
the ion acceleration can occur in a regime, when the laser 



pulse quickly sweeps significant part of electrons away, 
|lL ll2Lll4l |. At high intensity, when the radiation pressure 
of the laser field become dominant, the plasma is accel- 
erated so that almost all the energy of the laser pulse 
is transformed into the ion energy, 5]. Below we show 
that, according to our simulations, in the case of thin 
targets and optimal laser pulse duration, the ion maxi- 
mum energy scales as the square root of the laser pulse 
power, as suggested in 14]. With increasing intensity, 
the RPD regime comes into play and the ion maximum 
energy becomes proportional to the laser pulse energy. 

The laser-driven ion acceleration also depends on the 
structure of the laser pulse. "Parasitic" parts of the laser 
pulse such as a pre-pulse and an amplified spontaneous 
emission (ASE) can substantially change the structure of 
the solid target before the main part of the laser pulse 
arrives, e. g., a pre-plasma can be formed. The effect of 
"parasitic" parts of the laser pulse on the ion acceleration 
and dependence of the ion energy on the scale-length of 
the pre-formed plasma was dicsussed in Refs. (lfL.16>.17j |. 

We study the simplest model, where the "clean" gaus- 
sian p-polarized laser pulse, whithout "parasitic" parts, 
is incident at right angle on a thin slab of ideal collision- 
less plasma. The approximation of coUisionless plasma 
can be satisfactory when a femtosecond laser at inten- 
sity > 10^°W/cm^ interacts with plasma slab, which is 
a few laser wavelengths thick, even at solid density of 
plasma. The laser pulse is characterized by the wave- 
length A, the dimensionless amplitude a, corresponding 
to intensity / = a^h, h = 1.368x101^ W/cm2x(/im/A)^ 
the length (FWHM) L and focal spot diameter (FWHM) 
D. The target consits of two layers, the first layer is 
fully stripped aluminium and the second layer is a pro- 
ton coating. The first and second layers of the target are 
described, respectively, by the electron density and 
ne2, the thickness I and I2, the transverse (perpendicular 
to the laser pulse direction) size w and W2- 

We carried out two-dimensional (2D) PIC multi- 
parametric simulations using a technique described in 
|l7l |. In this technique a series of 2D PIC tasks with N 
different sets of the laser and target parameters is per- 
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formed simultaneously on the N processors of a multi- 
processor supercomputer, using the massively-parallel 
and fully vectorized code REMP, based on the PIC 
method and the "density decomposition" scheme [iSj. 
We analyze the dependence of the interaction outcomes 
such as reflection, transmission and absorption coeffi- 
cients, maximum ion and electron energy and kinetic 
temperature, beam emittance, energy spread, the accel- 
eration time and length, etc. on parameters of the laser 
pulse and the target: laser pulse intensity /, focal spot 
size D and duration L/c, target density and thickness 
I. In the simulations, the laser pulse initially propagates 
along the x-axis in the simulation box with x x y size 
251A X 136A; the transverse size of the target first layer 
is fixed: w = 80A, for the second layer it is W2 — D/2; 
the second layer is 0.06A thick, its density is such that 
the number of ions in the first layer in the longitudinal 
direction is ~ 10'^ times greater than the corresponding 
number of protons. At chosen conditions, protons behave 
almost like probe particles, thus we can analyze, almost 
independently, both the proton acceleration and the ac- 
celeration of heavier ions from the first layer. Simula- 
tions were performed on 720-processor HP AlphaServer 
at JAERI-Kansai. In each processor the grid consists 
of 4016 X 2176 cells and a number of quasi-particles is 
1.9 X 10^, amounting to 6.3 x 10^ of cells and 1.4 x 10^° 
of quasi-particles in total. 

Fig. n shows how the maximum proton energy de- 
pends on the target thickness and density and the laser 
intensity for fixed laser length L = lOA and focal spot 
size D — lOA. The proton energy was taken at the time 
when the proton acceleration dropped to 3% from max- 
imum; the corresponding acceleration length is approx- 
imately equal to the laser focal spot size. Frames (a)- 
(e) represent results of 720 separate computational tasks 
performed simultaneously in one run on 720-processor su- 
percomputer, each point with coordinates (l,ne) at fixed 
intensity / corresponds to particular one-processor task. 
The target thickness changes in the range Z/ A = 0.1 — 5 
(12 samples), the target density varies from rie = ricr to 
lOOricr (12 samples), where ricr = '"■/(^eA^) is the critical 
density, Ve — e'^/meC^ is the classical electron radius. The 
intensity varies from / = 10^° to lO^^W/cm^ x (^m/A)^ 
(5 samples) . Plasma slab with comparatively low density 
and large thickness can be considered as a model of the 
target modified by "parasitic" parts of the laser pulse. 

In fig. n the ion energy curves appears to be nearly 
parallel to the contours of the product a = nj, i. e. the 
number of electrons in the longitudinal direction. We 
call this parameter the critical depth of the plasma slab. 
We see that the energy dependence on two parameters 
I and rie degenerates to a dependence on only one pa- 
rameter a. Such a degeneracy is rather surprising if we 
inspect the results of different simulations with the same 
a. At small thickness / and large density Ug the laser 
pulse sweeps away a substantial part of electrons, the in- 



duced strong Coulomb potential of the first layer acceler- 
ates protons, in accordance with the scenario described 
in Refs. [EllJillJI- large thickness and small density 
the laser pulse penetrates throuhg the target and gener- 
ates a strong quasi-static magnetic field whose pressure 
causes charge separation which accelerates protons near 
the plasma-vacuum interface, similarly to the mechanism 
suggested in Refs. 0, 0|. The maximum energy gain 
corresponds to some optimum critical depth aopt- For 
decreasing a < Uopt, more and more laser pulse energy 
is transmitted through the plasma slab. For increasing 
c > fopt, the laser pulse reflection becomes more and 
more efficient. The energy absorption turned out to be 
optimal for the ion acceleration ai a — a^pt- Fig. ^f) 
shows the dependence of the energy gain on the laser in- 
tensity / and the critical depth a. Varying the ct, one 
can find the minimum intensity which gives the desired 
energy gain. In other words, at each intensity there is an 
optimal cr which gives the maximum energy gain. The op- 
timal critical depth almost linearly depends on the square 
root of the laser intensity, (7opt/"-crA « 3 + 0.4 x (I jl])^!'^ . 
This is similar to the criteria of relativistic transparency 
of a thin foil, derived in Refs. [20j . 

Fig. 121 shows the maximum proton energy vs the laser 
length and the target thickness for different laser in- 
tensities and focal spots at fixed plasma density rie = 
lOOricr. As above, the proton energy corresponds to the 
time when the acceleration damped to 3% from max- 
imum. The laser pulse length is chosen in the range 
L/A = 10 — 60 (6 samples), the target thickness is 
IjX = 0.125 — 5 (8 samples), the laser intensity sam- 
ples are the same as above, the laser focal spot diameter 
choices are DjX = 10,25,50. We see that the energy 
gain increases with decreasing target thickness and in- 
creasing laser pulse length. Columns (b) and (c) indicate 
that energy in the range 100-200 MeV is achievable with 
a "clean" petawatt laser with sub-picosecond duration 
when the solid density target is sufficiently thin, as was 
found in [13 . 

At intensity of the order of 5 x 10^-^W/cm^ x (^m/A)^, 
the new regime of acceleration comes into play, when the 
energy gain rapidly increases with the laser pulse length 
and, for optimal pair of laser duration and target thick- 
ness, relativistic protons can be obtained. This is the 
radiation pressure dominant regime of the ion accelera- 
tion, described in for much higher intensity. In this 
case the laser radiation pressure dominates in the interac- 
tion and the effective cross-section of the process (analog 
of the Thomson cross-section) becomes 2/nel, p. As 
shown in Fig. |21 the energy gain decreases if the laser 
duration is greater than the optimum. This is because 
we use Gaussian laser pulses, so a long pulse head have 
enough time to deteriorate the target. As indicated in 
Refs. 0, in the ideal case of the RPD regime, the energy 
gain is proportional to the laser pulse duration. 

Figs. 121 and 21 are compiled from the results presented 
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above. In Fig. O the maximum proton energy fmax is 
drawn for every instance of the laser pulse energy 
in the case I — X, = lOOnc,-. Points become arranged 
along 3 lines corresponding to 3 versions of the focal spot 
size; lines can be fitted by scaling fmax oc £^^^ up to 
^^max ^ 200 MeV and by fmax oc for higher energy 
gains when begins a transition to the RPD regime. 

Fig. 01 shows a correlation of the maximum proton 
energy with the laser power V for Ue = lOOncr, where 
points correspond to a thickness I close to optimal, i. e. 
Uel « CTopt. Narrowing the set of points by additional 
constraint that the laser pulse length is of the order of 
the focal spot size, L ~ Z), we obtain a strip aligned along 
the dependence £max « "P^/^. Once again we see a de- 
generation of a complex dependence of the ion maximum 
energy from many parameters to a dependence from a 
fewer number of parameters. This can be a manifesta- 
tion of a distinct regime. In this case the conditions are 
met for the regime described in Refs. P, 0, The 
square-root dependence on the laser power is consistent 
with the prediction made in Ref. Following Refs. 

P, 0, Q|, one can represent the irradiated spot as an 
uniformly charged oblate ellipsoid with size Ix D x D, 
I < D, and charge density +r]ene, where 7y is a portion of 
electrons swept away by the laser pulse. Assuming that 
the target critical depth is optimal, Ugl = OAaUcrX, one 
can obtain the maximum energy of the ion with charge 
Ze accelerated by the electric field of the ellipsoid (for 
I < D): £:max ~ riZ{P[PW])^/^ x 228 MeV, where V is 
the laser peak power. The condition that the optimal 
laser pulse length should be greater than or of the order 
of the focal spot size, L > D, turs out to be necessary 
to prevent return currents from affecting the ion acceler- 
ation (both the current, formed by returning electrons, 
and that, produced by electrons from a surrounding less 
irradiated and relatively cold plasma). 

In conclusion, our multi-parametric simulations show 
that in the process of the ion acceleration by an intense 
laser pulse from a double layer target, at given laser in- 
tensity / the dependence of the ion maximum energy 
from the target thickness I and density Tig is reduced to a 
dependence on the critical depth a = nj,. At given inten- 
sity the highest ion energy gain occurs at optimal critical 
depth of the target tJopt which is approximately propor- 
tional to the square root of intensity, (Jopt oc 1^1'^ . If the 
target critical depth is less than the optimal value, the 
laser pulse is more transmitted rather than absorbed; if it 
is greater, the laser pulse reflection increases making the 
interaction less efficient for the ion acceleration. Simula- 
tions indicate that in the case of thin targets and optimal 
laser pulse duration, the ion maximum energy scales as 
the square root of the laser pulse power, in consistence 
with the theory presented in Ref. Therefore, the 

laser-driven ion acceleration up to 200-300 MeV, which 
is necessary for hadron therapy, can be obtained with a 
petawatt laser with sub-picosecond duration. A transi- 



tion to the radiation pressure dominated regime of the 
ion acceleration is revealed. With increasing intensity, 
when the RPD regime comes into play, the ion maximum 
energy becomes proportional to the laser pulse energy. 
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FIG. 1: Maximum proton energy (contours) vs. target tliicltness and density (log-log scale) for different laser intensities, 
(a)-(e), for L — lOA, D = lOA. The laser pulse reflection (di), absorption (d2) and transmission (da) coefficients, related 
to frame (d). Dashed line in (d), (di-a) - for the best Uel, corresponding to maximum energy gain, (f): maximum proton 
energy vs laser intensity and target critical depth a = Uel (log-log scale). Dashed line for optimal aopt oc I^^^. Intensity unit 
W/cm^ X i/iia/Xf. 




FIG. 2: Maximum proton energy (contours) vs. laser length and target thickness (log-log scale) for different laser intensities 
(columns (a)-(e)) and different laser focal spots (raws 1-3), for rie — lOOricr- 
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FIG. 3: Maximum proton energy vs. laser pulse energy for FIG. 4: Maximum proton energy vs. laser power for optimal 
I = X, ris = lOOiicr- Dashed lines exemplify possible scalings. plasma slab thickness and rie = lOOricr- 



